Light cones in relativity: Real, complex and virtual, with 

applications 



T.M. Adamo 
The Mathematical Institute 
University of Oxford 
24-29 St Giles' 
Oxford, 0X1 3LB, U.K. 
adamoOmaths . ox . ac . uk 



E.T. Newman 
Department of Physics & Astronomy 
University of Pittsburgh 
3941 O'Hara St. 
Pittsburgh, PA 15260, U.S.A. 
newmanOpitt . edu 



January 6, 2011 



Abstract 

We study geometric structures associated with shear-free null geodesic congruences in 
Minkowski space-time and asymptotically shear-free null geodesic congruences in asymp- 
totically flat space-times. We show how in both the flat and asymptotically flat settings, 
complexified future null infinity, 3^7, acts as a "holographic screen," interpolating between 
two dual descriptions of the null geodesic congruence. One description constructs a com- 
plex null geodesic congruence in a complex space-time whose source is a complex world- 
line; a virtual source as viewed from the holographic screen. This complex null geodesic 
congruence intersects the real asymptotic boundary when its source lies on a particular 
open-string type structure in the complex space-time. The other description constructs a 
real, twisting, shear-free or asymptotically shear-free null geodesic congruence in the real 
space-time, whose source (at least in Minkowski space) is in general a closed-string struc- 
ture: the caustic set of the congruence. Finally we show that virtually all of the interior 
space-time physical quantities that are identified at null infinty 3^, (center of mass, spin, 
angular momentum, linear momentum, force) are given kinematic meaning and dynamical 
descriptions in terms of the complex world-line. 
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1 Introduction 

In this paper, we describe some interesting structures based in classical special and gen- 
eral relativity which bear some resemblance to dualities known as "holographic" dualities 
which have emerged elsewhere in theoretical physics over the past decades (c.f , [1, 2, 3]). 
Though these holographic dualities usually involve the use of highly non-classical machin- 
ery such as supersymmetry or string theory, most famously in the case of the AdS/CFT 
correspondence [4, 5, 6], we emphasize that our discussion here will use no such tools; 
we work entirely in the context of classical four-dimensional Lorentzian space-time. The 
structures we are interested in emerge naturally from the study of light cone foliations (and 
their generalization to asymptotically shear-free null geodesic congruences) in space-time, 
and have simply been overlooked in prior research. Although it would certainly be pre- 
sumptuous for us to suggest that our work here has any true connection with holographic 
duality as it is known to most theoretical physicists, we do find a holographic screen, open 
and closed classical strings and other suggestive objects; all of which can be given real 
physical meaning in four-dimensional space-time. 
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Specifically, it is the purpose of this note to first explore the properties of ordinary "run 
of the mill" light-cones and then turn to their generalization via complex and virtual light 
cones in four-dimensional Lorentzian space-times. More precisely, we study the properties 
of light-cones and their complex generalizations in both Minkowski space and in asymp- 
totically flat (vacuum and Einstein-Maxwell) space-times in the neighborhood of future 
null infinity. This is followed by a discussion of physical applications of these ideas and 
constructs. These complex light cones are first applied to the structure of real Maxwell 
fields in real Minkowski space. The complex cones in flat space-time are then generalized 
and applied to study equations of motion in general relativity. Along the way we point 
out a pretty duality between the complex light-cones and real shear-free, but twisting, null 
geodesic congruences. 

The first issue raised comes from the simple question: In Minkowski space, avoiding 
or ignoring their apex, what are the geometric properties that a set of null geodesies 
must have in order to be a light cone? How can they be determined to be light cones 
even far from their apex? The answer is simple: First of all, the family of relevant null 
geodesies (the light-cone generators) must be (null) surface forming; they must lie on a 
null surface and thus have vanishing "twist." As a result, these surfaces must be foliated 
by null geodesies whose tangent vectors are determined by the gradient of the surface. 
Second, they must have vanishing shear and non-vanishing divergence. (The plane null 
surfaces can be thought of as light-cones but with their apex at infinity; we ignore this 
case.) From this requirement, light cones posses the topology of x K. For us, the most 
relevant feature is their vanishing shear. Even far from their apex (i.e., even at future null 
infinity, 3"*") on the portion of the surface, if the shear vanishes then it has an apex 
and the surface is a light-cone. This case can be generalized from an individual light-cone 
to a family of light-cones: if Minkowski space is (partially) foliated by a null geodesic 
congruence (NGC) , can we tell at that the geodesies all focus to a time- like world-line 
in the interior? The answer again is simple: If the congruence has vanishing twist and 
shear at 3+ and non-vanishing divergence, and furthermore has no members lying tangent 
to 3"*" itself (the regularity condition) , then there is a one-parameter family of light-cones 
and a time-like world-line to which the NGC converges. 

The main goal of this work is to investigate and analyze how this asymptotic description 
of light cones can be generalized, and what applications to physics it might have. The 
generalization will be in two distinctly different but related directions. 

First of all, in the context of Minkowski space, we define and describe complex light 
cones. They will be determined solely from the properties of specific sets of null directions 
at complexified null infinity (3^ ); analytic continuation of Penrose's future null infinity, 
3+. These complex null directions, normal to specific slices of 3^, define - by following 
them backwards in time - complex null geodesies (and complex light-cones) which converge 
to points in complex Minkowski space [7, 8]. In general there will be a subset of points 
in complex Minkowski space where one (or more) of its light-cone generators intersects 
real Minkowski space at real 3+. If instead of the complex light-cone of a single point in 
complex Minkowski space we take a complex analytic "time-like" (to be defined) world-line 
parametrized by the complex parameter r, we would have a two-real-dimensional set of 
complex points (from the real and imaginary parts of r) and their light cones. We show 
that for any fixed value of the real part there is a one-dimensional set of points such that 
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the envelope formed by their individual light cones intersects real 3+ on an slice. As 
the real parameter changes we obtain a one parameter family of real slicings of 3^ . 

At this point, an interesting duality emerges. On one hand, if we start from each of 
these slices and move backwards into the complex space along the complex null directions, 
these trajectories converge to an imaginary line segment in the complex space. On the 
other hand, there is a dual method (described later) for following null geodesies from the 
slices back into the real space-time; this yields a real shear-free, but twisting, null geodesic 
congruence. It is precisely this twist which links the two pictures: the "distance" of the 
complex world-line from the real Minkowski space-time in the first picture is a measure of 
the twist of the real congruence in the latter picture. The caustic set of the real (dual) 
congruence is (in general) a closed curve moving in real time; something analogous to a 
classical closed string [8]. 

The extension of these ideas to asymptotically flat Einstein space-times initially seems 
to be impossible. Standard light cones from any given space-time point will undergo such 
distortions from the curvature of the space-time itself that little or no memory of their 
origin will remain when they arrive at . Nevertheless, we can consider the possibility 
of using the procedure that was succesful in the Minkowski space-time case by asking for 
null geodesic congruences in the neighborhood of that are shear-free and non-twisting. 
In the general asymptoically flat case, shear-free null geodesic congruences do not exist 
- but there are always null geodesic congruences that are asymptotically shear-free in 
the neighborhood of J"*". Unfortunately, to use this idea effectively again entails the 
analytic extension of the space-time a small distance into the complex. Working on the 
complexification of 3'^ (i.e., on 3^), there is a construction of complex slices or "cuts" 
whose complex null normals can be used to determine asymptotically shear-free and twist- 
free complex null geodesies [9, 10, 7]. In fact one can construct a four-complex dimensional 
family of such complex cuts which define a four-complex dimensional manifold frequently 
referred to as H-space [11, 12, 13].^ The immediately relevant feature for us is that these 
complex null geodesies from each complex "cut" converge or focus to a point in H-space 
[14] . It will be shown later that real structures associated with H-space can be found, and 
that real physics can be interpreted as taking place in H-space [15, 7]. The "H-space can 
thus be viewed as the virtual image space seen by looking-backwards along complex null 
directions from a sphere of points on 3^. It is this property that could allow us to refer 
to 3j as a holographic screen. 

The prior discussion of complex Minkowski space (which is a special case of "H-space) 
can be extended to "H-space. There is a subset of points in 7^-space where one (or more) 
of the light cone generators (null geodesies) coming from a complex point intersects real 
asymptotically flat space-time at real If instead of the complex light cone of a sin- 
gle point in complex "H-space, we take a complex "time-like" (to be defined) world-line 
parametrized by the complex parameter t — s -\- iX, we would have a two-real-dimensional 
set of complex points (from the real and imaginary parts of t) and their associated complex 
light cones. For any fixed value of the real part, s, there is a one-dimensional set of points 
(a finite interval parametrized by A) such that the envelope formed by their individual 

^Though it is not used here, we mention that "H-space is endowed with a complex Ricci-flat metric with 
self-dual Weyl tensor. 
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light-cones intersects real 11+ on a S"^ slice. As the real parameter s changes we obtain in 
the "H-space a ribbon (the finite interval moving in "s-time" ) which could be called a classi- 
cal open string; from the null cones of points on this ribbon, we get a one parameter family 
of real slicings on 3+. All the information about the ribbon is encoded (holographically) 
in the one parameter family of real 3+ slicings and a null direction field on In other 
words, there is a duality between the coded information on 3"*" and H-space information. 
A further related duality is that a given complex analytic world-line in "H-space (via its 
associated ribbon) yields in the physical space-time (via a complex conjugate action) a real 
twisting but asymptotically shear-free null geodesic congruence in the real space-time. 

The question of where this beautiful mathematical structure makes contact with phys- 
ical issues does have a simple answer - the details, however, are rather complicated. 

The simple answer is that an (analytic) asymptotically flat Maxwell field in Minkowski 
space with non-vanishing total charge generates, in the complex Minkowski space, a unique 
complex analytic world-line: the complex center of charge world-line, where the real part 
describes the standard center of charge while the ribbon thickness encodes magnetic dipole 
information [16, 17]. For the case of asymptotically flat space-times there are two situa- 
tions: the vacuum asymptotically flat and the Einstein-Maxwell asymptotically flat space- 
times. For the vacuum case there is a unique complex H-space world-line that contains: 
from the real part, the equations of motion for the physical center of mass and, from the 
ribbon thickness, the spin angular momentum, with both interpretations arising from the 
"view" at infinity, 3''" [15, 7]. These are (loosely) analogues of measuring the total charge 
at infinity via Gauss's law or observing the Bondi energy-momentum vector at infinity. 

In Section 2, the preliminaries, we introduce our notation and results from earlier in- 
vestigations that will be needed here. Specifically we first discuss conventions and notation 
followed by a description of flat-space null geodesic congruences. The section ends with a 
brief summary of properties of asymptotically flat spaces and their null geodesic congru- 
ences. Section 3 deals with real space-time structures that are associated with the complex 
world-lines; first in complex Minkowski space and then in H-space. In Section 4 we apply 
the ideas associated with the complex world-lines to real physical ideas. In particular we 
show that a real asymptotically flat Maxwell field (with non- vanishing charge) determines 
a complex world-line (the complex center of charge) that carries information about both 
the electric and magnetic dipole moments. This construction is then generalized to asymp- 
totically flat space-times, where the complex mass dipole moment (the real mass dipole 
moment plus "i" times the angular momentum) determines an H-space world-line. The 
Bianchi identities then yield kinematic definitions, equations of motion, angular momen- 
tum and conservation laws. All take place in "H-space which we interpret as a virtual image 
space. This information is coded into the real space-time by functions on real The 
results are very reminicent of ordinary Newtonian dynamical laws of motion. Though 
partially a summary of results presented elsewhere in the literature (e.g., [7, 15, 17]), our 
presentation includes several simplifications and alterations. In Section 5 we summarize 
the earlier discussion and speculate on what meaning and possible future use there might 
be to the observations made here. Appendix A provides some background on tensorial 
spin-s spherical harmonics, which are used throughout this work. 

We again stress that the strange results described here lie wholey in standard four- 
dimensional classical physics. There is no need - other than assumed analyticity - to 
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rely on drastic modifications of space-time properties such as supersymmetry or higher 
dimensions. The resuhs are here to be seen and perhaps understood. It would have been 
a cruel god to have layed down such a pretty scheme and not have it mean something 
deep. 

2 Foundations: Null Geodesic Congruences and 

Asymptotic Flatness 

In this section we summarize several of the basic ideas and tools which are needed in our 
later discussions. The explanations are rather concise and extensive proofs are omitted. 
In large part, much of what is covered in this section should be familiar to many or even 
most workers in general relativity. 

2.1 Conventions and Notation 

The arena for most of our discussion is the neighborhood of the "far (infinite) null future" 
of our space-time, (intuitively the end-points of future directed null geodesies) for both 
Minkowski space and asymptotically flat space-times. This region, first defined and studied 
by Roger Penrose and referred to as future null-infinity, 3"*" , is constructed by the rescaling 
of the space-time metric by a conformal factor which appoaches zero asymptotically, the 
zero value defining 3+ [18, 19, 20]. This process leads to the (future-null) boundary 
being a null hypersurface for the conformally rescaled metric with topology S'^ x R. An 
easy visualization of the boundary 3"*" is as the past light cone of the point future 
time-like infinity. A natural coordinatization of 3+ and its neighborhood is via a Bondi 
coordinate system: (u, r, C, C,). In this system, u, the Bondi time, labels the null surfaces 
of the space-time that intersect 3'^; r is the affine parameter along the null geodesies of the 
constant u surfaces; and C = e"^ 001(61/2) is the complex stereographic angle that labels 
the null geodesies of 3+ , the portion of 3"*" [9] . 

In Minkowski space, the Bondi coordinates (u, C,, () of 3+ can be constructed from the 
intersection of the future null cones of the time-like world-line at the Minkowski space 
spatial origin, i.e., from the line, = (t, 0, 0, 0). The cone has the form 



= UretS^ +rl- iCX), 

( = e'"^ cot (61/2), Uret^t-r^V2u 

- + C + C"-*(C -<),-! + CC) 

a/2 

= -^(I,sin0cos0, sin^sint/), cos6') 

with (C, C) labeling the sphere of null directions at the origin and Uret the retarded time. 
3^ is the limit as r tends to infinity. 
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Remark 2.1 Note that Uret,u,t (and the variable t introduced later) all have the dimen- 
sions of length. In the Section 4, the velocity of light, c, will be explicitly introduced via the 
replacement (uret,u,t, t) ^ {cUret,cu,ct^CT) so that Uret,u,t and r have the dimensions 
of time. 

Remark 2.2 We note that the round sphere metric, ds^ = dO^ + s\v? Odip^ becomes in 
stereographic coordinates ds^ = 4P^^d(^dC, with P — 1 + (^^. 

To reach 3+, we simply let r — oo, so that 3+ has coordinates (u, C, C)- The choice 
of a Bondi coordinate system is not unique, there being a variety of Bondi coordinate 
systems to choose from. The coordinate transformations between any two are know as 
Bondi-Metzner-Sachs (BMS) transformations or as the BMS group (c.f., [21, 22]). 

Our assumption of the analyticity of the space-time then allows for the complexification 
of J"*". For this complexification (i.e., extension to 3^), we allow u to take on complex values 
clos^ to the real and free ( from being the complex conjugatae of (. It is then denoted 
by C ~ C- (Often we take this as implicitly understood and just use C-) 

Associated with the Bondi coordinates is a (Bondi) null tetrad system, (Z", n", m", m") 
(c.f., [9, 23]): 



The first tetrad vector Z° is the tangent to the geodesies of the constant u null 
surfaces given by 



r = —=g-'^V,U, (2.1) 

rVa/' = 0, (2.2) 

The second null vector is tangent to the null geodesies lying on normalized so that: 

Ln" =1. (2.4) 

The remaining vector m° and its complex conjugate are tangent to the S"^ slices of constant 
u. 

An important construct is the family of past light-cones from each point of 3+ (or 3(t). 
Each past cone is determined by a sphere's worth of null directions at 3+ with each null 
direction labeled by the associated sphere coordinate. These coordinates are chosen as the 
complex stereographic coordinates and are denoted by the complex conjugate pair (L, L). 
An arbitrary field of null directions on (and consequently an arbitrary null geodesic 
congruence that intersects 3+) can then be described by the function L = L{u, () or its 
analytic extension to 3^^. 
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Often we will use a very specific form of a null tetrad given in Minkowski coordinates 
and parametrized by the points on the sphere in stereographic coordinates (C, C) and 
denoted by the over-hat: 

t ^ ^ (1 + cc, c + c -.(c -<-),-! + cc) ^ (f , inf ) , 



^ (o,i-c',-*(i + c'), 2C) = (o,-rW), 



2(1 + CC) 

V2 



(o,i-c2, z(i + c'),2C),= (o,-yi(-^^ 



2(1 + CC) 

As (C, (^) move over the complex plane, and n° range over the light cone. The 
spin-s harmonics [24], Yj^f-j, (C, (^), which are frequently used, are described in Appendix 
A. 

2.2 Flat space Null Geodesic Congruences 

In Minkowski space M, the future light-cones from an arbitrary time-like world-line = 
£,°'{s) can be described by the Null Geodesic Congruence (NGC) 

x'^=r(s)+rr(C,C) (2.5) 

with r the affine parameter on each of the light-cone generators. This construct is easily 
generalized to complex Minkowski space Mc, where light-cones from 5°(r) (now an arbi- 
trary complex analytic world-line with complex affine parameter r) and its corresponding 
complex NGC is ^ 

j'^ = C{r)+rr{Q,0, (2.6) 
where r is now complex and (C, C) ^re independent of each other. 

The Sachs complex optical parameters for an arbitrary NGC (real or complex) are the 
complex divergence and shear of the congruence [25, 23], 

P = ^(-V„r + i curU'^), (2.7) 
a = V^al^m^m^ (2.8) 

where 



curl r= ^(V[,?f,]V-/''). 
These satisfy the flat-space optical equations: 

Dp ^ + aa (2.9) 

Da = {p + p)a, 

D = rA.A, 

dx'^ dr ' 
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with r the affine parameter along the geodesies. The optical parameters for the above 
light cone congruence can be calculated directly from Eq. (2.5) yielding 

P = -r-\ (2.10) 
a ^ Q. 

By reversing the statement and assuming a NGC with vanishing shear and real diver- 
gence, the optical equations become 

Dp = p\ 

The integral (i.e., p — ~r~^) is the same as Eq.(2.10), thus showing that a NGC with real 
divergence and vanishing shear is the light-cone congruence of a (real) time-like world-line. 

An arbitrary NGC in Minkowski space can be described by the three parameter, 
(u, C, C)i family of null geodesies 

= uil" + n") - Lm" - Lm + (r - ro)t (2.11) 

where r is the afhne parameter, tq = ro(w, CjC) is the arbitrary origin for the afhne pa- 
rameter and L = L{u,(^,(), the determining function of the congruence, is an arbitrary 
complex spin-weight one function of the parameters. The three parameters (w, C, C) are 
the Bondi coordinates of the intersection points of the null geodesies with 3+. The op- 
tical parameters are determined by L which is the stereographic angle field on 3+ that 
determines the directions of the null geodesies. 

The condition for a NGC with vanishing shear is that the function L must satisfy the 
non-linear partial differential equation [26] 

dL + LL = 0, (2.12) 

where L ~ L,„, and 3 is the spin- weighted covariant derivative on the 2-sphere (see 
Appendix A for details) [27]. This can be integrated by introducing an auxiliary complex 
variable r = T{u, (, C), related to L by the so-called CR equation (related to the existence 
of a CR structure on 3+ [28]), 

gr + LT = 0, (2.13) 

and then using its inversion 

u = G(r, C,C),_ _ (2.14) 
T - r(G(T,C,C),C,C) = T. (2.15) 

After a process of implicit differentiation (c.f., [29, 10, 7]), Eq.(2.12) becomes 

g^G^O, (2.16) 

with 

the subscript r indicates that the differentiation is at t held constant. From this it follows 
that the regular solutions to Eq.(2.12) can be given implicity in terms of G as: 
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u = G{TX,0 = Cir)la^T=JiuX,0, (2.17) 

Several remarks must be made here: 

• are four arbitrary complex analytic functions of the complex parameter r which 
can be interpreted as determining a complex world-line in complex Minkowski space. 

• r can be regauged by the analytic function t* — F{t). Often it is useful to chose 

• Since r is complex we must allow u to take complex values which requires the com- 
plexification of 3+ , denoted 3^ . 

• When the ^°(t) are real functions of a real variable, s, Eq.(2.11) reduces to Eq.(2.5) 
(i.e. to the real world-line light cone congruence). 

• For a complex set ^°(t), the NGC, Eq.(2.11), is a real shear-free NGC but with a 
non-vanishing twist. The caustic set is in general a closed curve moving in time. 

An important observation that plays a major role for us is the following: From the same 
L, two different "conjugate" versions can be constructed. The first is obviously the com- 
plex conjugate given by L = ^ {T)Tha while the second, referred to as the holomorphic 
conjugate, is given by L = SG = ^°(r)TOa. Using L in Eq.(2.11) instead of L, we obtain 
another shear free NGC but now it is the complex congruence, given earlier by Eq.(2.6), 

^" = r(T) + rr(C,C). 

In other words, the cut function describes a family of null cones with apex on the 
complex line, z° = £,'^{t). We now have on 3^ and 3^ two different tetrad systems 
(obtained by null rotations from the Bondi tetrad) coming from L and L, namely 

p ^ l*^ ^kra'' --m" + 0{r~'^), (2.18) 
r r 

m*" = m"- -n", (2.19) 



r 



(2.20) 



and 



r ^ l^" =r--m''--m'' + 0{r-^), (2.21) 



r r 



-n°, (2.22) 



r 



(2.23) 



The null geodesic congruence determined by as mentioned earlier, is a real shear- free 
congruence with twist while the congruence determined by 1^°' is complex, shear-free, twist 
free conguence and focuses on the complex curve f (t). 
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Though the complex nuh geodesies with apex on ^"(t) spend most of their "time" in 
the complex Minkowski space some do reach real Minkowski space and in particular some 
reach the real . It turns out that the complex world-line and their associated light-cones 
have real structures. They are discussed in Section 3. 



2.3 Asymptotic Flatness 

At a first glance it would appear as if it were not possible to duplicate the Minkowski space 
discussion of light-cone NGCs in asymptotically flat space-times. Aside for a few special 
cases (the algebraically special metrics) there are no Einstein space-times with shear-free 
NGCs. The family of future directed null geodesies originating at a fixed space-time point 
traversing regions of curvature will, in general, be distorted and develop shear. Surprisingly 
it nevertheless is possible to duplicate virtually all the light-cone NGC results of flat- 
spacetime for the general case of asymptotically flat space-times by looking not for shear- 
free NGCs but instead asymptotically shear-free conguences. In fact such congruences are 
determined by a complex analytic curve in an auxiliary four-complex dimensional space, 
referred to as "H-space. 

Before describing these congruences we first review some relevant features of asympot- 
ically flat space-times. Details, derivations and proofs are largely omitted since they are 
easily found in the literature [25, 23, 9, 7]. 

We begin by pointing out that with Bondi coordinates and tetrad the two optical 
parameters, the complex divergence and shear are given by 

p = p=-i + ^ + 0(,-5), (2.24) 
a = ^ + 0{r'% 

with = fT"(w, C, C)j the asymptotic Bondi shear of the NCG with the Bondi tangent 
vector, i.e., l"". The cr°, which is the free data determining the gravitational radiation, plays 
a major role in our discussion. Considering a new NGC with tangent vector /*" defined 
at 3+ by the null rotation 



Z*'^ = r - -m" - -m'' + 0{r-'^), (2.25) 
r r 

m*° = m''~-n'\ (2.26) 
r 

= (2.27) 

with arbitrary L = L(w, C, C), one finds that the asymptotic shear of the new conguence is 
given by a version of the Sachs theorem [26] : 

a°* =BL + LL- a°. (2.28) 

The condition for the new congruences to be asymptotically shear-free ((T°* — 0) is thus 
that L satisfy 

dL + LL^a°, (2.29) 
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which is the extension of the flat-space Eq.(2.f2). 

As in the Minkowski space case, Eq.(2.12), this can also be integrated by introducing 
the same auxiliary complex variable r = T(u, C, related to L by the CR equation, 

dT + Lf = 0, (2.30) 

and then using its inversion 

u - G(r, C,C),_ _ (2.31) 

T = r(G(T,C,C),C,C) = T. (2.32) 

Note that as in the flat case, r is complex and we must allow the complexification of 
u and let C ~ C- For each value of r we obtain a complex "cut" of 3^. 

Again after manipulating several implicit derivatives, Eq.(2.29) and (2.30) become 

dlG = f7°(G,C,C), (2.33) 
L(u,C,C) - ^rG\^=TiuXX) (2-34) 

with again the subscript t indicating that the derivatives are at r held constant. Eq.(2.33), 
the "good cut equation," has been shown to depend on four complex parameters, (the 
H-space coordinates), so that we can write, u = X{z°-, C, (), where we distinguish X from 
G by its explicit dependence on the four solution parameters. By the coordinate freedom, 

the first four spherical harmonic coefficients can be chosen as the z"', (coordinate conditions 
on H-space) so that we have 

M = X(z^ C, C) = ^''^a (C, C) + Hi^2 (z% C, C) , (2.35) 

where Hi>2 are spherical harmonic contributions with I > 2. Finally, by taking an arbitrary 
world-hne in the "H-space, z° = C'^(r), we find the general regular solution to Eq.(2.29) is 
given implicitly by 

L(u,C,C) = ^rGl^^^^cro (2.36) 
G(r, CO ^ X(r(T),C,C) (2.37) 
u = G(r, C,C) (2.38) 

= r(T)/;(c,c)+i?;^2(r(T),c,c) = r(T)fa(c,c)+%2(T,c,c). 

As in the flat space case, the regular asymptotically shear-free NGCs are determined 
by the arbitrary choice of a complex world-line in an auxiliary complex space, ?^-space. 

In complete analogy with the complex Minkowski space case, it turns out that if we use 
the complex NGC determined by the stereographic angle field, (2.36), and the associated 
holomorphic field, L(u, C,C) = ^rG\^^j.^^^^ -^y (not L(u, C,C) = Si-G|:^^j^(„ q) as initial 
directions at 3^ 

l*^ =r - Lm''^ -m" + Oir-^), (2.39) 
r r 



12 



the complex geodesies converge on the 'H-space world-Hne z°- — C(t) [14]. By pointing 
into complex null directions from the complexified 3^7 , we have complex virtual cones and 
complex virtual world-lines. There will always be points on the complex world-line whose 
null cones partially intersect real 3+. We will see later that unique world- lines can be 
determined so that real meaning or significance can be given to them, as complex centers 
of charge and complex centers of mass - which include both asymptotic magnetic dipoles 
and angular momentum. 

The basic idea that will be pursued later (in Sections 3 and 4) is to identify certain 
terms in the asymptotic behavior of the Maxwell and Weyl tensor tetrad components 
with physical quantities and then see how they change when they are computed with 
the (rotated) complex null directions pointing towards a complex world-line, Eq.(2.39). 
By choosing the world-line appropriately, so that these quantities vanish, we identify the 
virtual complex centers of charge and mass. 

As an interim step, we need the behavior of both the tetrad components of the Weyl 
and Maxwell tensors: 

^0 = -Cabcd^'m^rm^ (2.40) 
^1 = -Cabcdl''n'l'ni'', 

4^2 = -^{Cabcdl''n''l-n''-Cabcdl''n''m'rff'), 

I /" i ja b c — d 

V3 = Cabcdl n n m , 
tpi = Cabcdn°'rrfri''rff'-, 

and 

00 = Fabl^m^, (2.41) 

01 = iF,b(rn'' + m"m^), 

Integrating the both the Weyl tensor and Maxwell spin-coefhcient equations leads to 
the peeling behavior: 



■00 


= <r- 


-5 _ 


hO(r- 


-% 


V'l 


= 0°r- 


-4 _ 


hO(r- 


-% 


■02 


= 02°^^ 


-3 _ 


f-0(r- 


-% 






-2 ^ 


hO(r- 


-% 


■04 




-1 


^0{r- 





(2.42) 



and 



+ 0(r-4), (2.43) 



0°r-2 + 0(-3N 



r 



62 = 0^r-i+O(r-2^ 
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With the coefficients satisfying the asymptotic Bianchi Identities 



^§ = 


- 




^? = 








-90? - 





(2.44) 
(2.45) 
(2.46) 



where 



and asymptotic Maxwell equations 



V'4° - (2.47) 



• 



V'3 = 9a", (2.48) 



0? - -90^, 



(2.49) 

(T>^'. (2.50) 



All the coefficients, which are functions on 3+ , i.e., functions of (m,C, i^), have 
physical meaning, e.g., multipole moments, etc. The four quantities {V'? , '02; 0oi 0i} 
the most importance to us, due to the following properties: 

• The I = harmonic of is the Coulomb charge q. It is assumed to be nonvanishing 
whenever a Maxwell field is being considered. 

• The / = 1 harmonic of 0q is the complex electromagnetic dipole moment: -D|;^jv,/ = 

• The Z = 0, 1 harmonics of -i/j^, slightly modified by the shear cr", is the Bondi energy- 
momentum four- vector. 

• The I = 1 harmonic of also slightly modified by the shear (T°, encodes the center 
of mass dipole and angular momentum: -Dcjg^Q.y) = D^^^^^ + icr^J\ 

By using the tetrad transformation generated by Eq.(2.39) (see Eq.(2.21)), one find 
the transformation law of the leading terms of the Weyl and Maxwell tensors: 

- iL^i;f + L^i'l (2.51) 
i^^?, (2.52) 

(2.53) 
(2.54) 
(2.55) 

C - 0°-2L0; + L2 0O, (2.56) 
? - ^02, (2.57) 

(2.58) 







4LVf + 6L^^/4" 




= V?- 


ZL^l + 3X^3 - 














^? 







Later, setting to zero the / = 1 parts of ■01*' and 00^, we can determine two different 
world-lines (when a Maxwell field is present) that can be referred to respectively as the 
complex centers of mass and charge. 
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3 Real Structures from the Complex World-line 



Our task in this section is to find the real structures that are lying in the complex world- 
lines and their complex light-cones. 

3.1 Flat-space Real Structure 

We first examine the case of flat-space-tinie with a complex Minkowski world-line, z° = 
^"(t) and associated light-cone cut of 3^, 

u = e{r)la{C,0- (3.1) 

To answer our question: what values of r allow real values of u, we first write r = 
s + iX, (s, A real) and decompose the right-hand side of Eq.(3.1) into its real and imaginary 
parts and set the imaginary part to zero [8], 

U = l{ei-^ + ^X)la{CX)+T{s-tX)laiCX)) (3.2) 

= [e{s + ^X)-T{s-^X)]la{CX)■ (3.3) 
Considering Eq.(3.3) as an implicit equation defining 

A-A(s,C,C) (3.4) 
we have that the allowed values of t are given by 

r = s + zA(s,C,C)- (3.5) 
The real values of u are thus given by the one-parameter (s) family of slicings 

u^eis + ^A{s,C,0)la{CX)■ (3.6) 

Assuming small values for the imaginary part of ^'^(t) — ^^(r) (r), (C^ ('''), ^i{t)) 
both real analytic functions) and hence small A(s, C,, (), it has been shown that A(s, (), (for 
fixed value of s), is a bounded smooth function on the (C, C) sphere, with maximum and 
minimum values, Amax = A(s, Cmax, Cmax) and Amin = A(s, Cmin_2_Cmin)- Furtherniore on the 
sphere, there are a circle's (S^) worth of curves between (Cmin, Cmin) and (Cmax, Cmax) such 
that A(s, C, C) is a monotonically increasing function on each curve. Hence there will be 
a family of circles on the (C, C)-sphere where value of A is a constant, ranging between 

Amax and Amin. 

Summarizing, we have the result that in the complex r-plane there is a ribbon or 
strip given by all values of s and line segment parameterized by A between Amin and 
Amax such that the complex light-cones from each of the associated points, ^"(s -I- iX), all 
have some null geodesies that intersect real . More specifically, for each allowed value 
of T = s -|- zA there will be a circle's worth of complex null geodesies leaving the point 
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^"(s + iX) reaching real 3+. It is the union of these nuh geodesies, corresponding to the 
circles on the (C, C)-sphere from the line segment, that produces the real family of cuts, 
Eq.(3.6). 

The real structure associated with a complex world-line is then the one-pameter family 
of slices (cuts) Eq.(3.6) and angle field L{u, (, Q on each point of the cuts. 

The dual point of view, as previously mentioned, is to start with the same L as used 
earlier: 



u = e{r)ia^r = T{u,(:,0, _ (3.7) 
L = 5,G = r(T)ma=r(n",C,C))ma. 

which was used with the holomorphic L, 

but now, instead, use the complex conjugate of L: 

L = drG = 1'\T)rha 

for the null directions pointing inward. In this case one obtains again a real shear-free 
NGC but now with twist S(u, C) which comes from the complex divergence, 

P = (3.8) 

r + i2j 

2iE = BL + L(Ly -dL-LL. (3.9) 
= iCir) -T if)) {na- la). 

As was claimed earlier, the twist is proportional to the imaginary part of the complex 
world- line and consequently we have the real structure coming from two (dual) places. 

3.2 Asymptotically Flat-Space Real Structure 

The extension of the above argument to the case of asymptotically flat space-times is 
relatively simple. Again assuming that the Bondi shear is sufficiently small and the 
■H-space complex world-line is not too far from the "real," the solution to the good-cut 
equation (2.33), 

u = r(r)fa(C,C) +i^;^2(r,C,C) = G(r, CO, (3-10) 
with T = s + iX, is decomposed into real and imaginary parts, 

G(t,C,C) = ^{G{s + iXX,0 + G{s-iXXX}) + l{G{s + iXXX)-G{s-iXXX))- (3.11) 

Setting the imaginary part to zero and solving for A we obtain an expression of the 
form, 

A = A(s,C,C)- 



16 



As in the flat case, for fixed s = sq, A has values on a line segment bounded between some 
Aniin and Aniax • The allowed values of t are again on a ribbon in the r-plane, all values 
of s and values on the A-line segments. 

Each level curve of the function A — A(so, Ci C) ^constant on the (C, C)-sphere, (closed 
curves or isolated points) , determines a specific subset of the null directions and associated 
null geodesies on the light-cone of the complex point ^"(sq + *A(so, C, 0) that intersect the 
real 3+. These geodesies will be referred to as "real" geodesies. As A moves over all allowed 
values of its segment, we obtain the set of H-space points, ^°(so + zA(so,C)C)) aiid their 
collection of "real" geodesies. From Eq.(3.11), these "real" geodesies intersects 3^ on the 
cut 

u =G(so+^A(so,C,C),C,C)• 
As s varies we obtain a one-parameter family of cuts. If these cuts do not intersect 
with each other we say that the complex world-line £,"'{t) is by definition a "time-like" 
line. This occurs when the time component of the real part of the complex velocity vector, 
v°'{t) — d^''(T)/dT, is sufficiently large. 

3.3 Summary of Real Structures 

To put the ideas of this section into perspective we collect the claims. 

• In Minkowski space, the future directed light-cones emanating from a real time-like 
world-line, a:" — ^°'{s), intersect future null infinity, 3+, on a one-parameter family 
of spherical non-intersecting cuts. 

• The complex light-cones emanating from a time-like complex analytic curve in com- 
plex Minkowski space, z° — ^"(t) parametrized by the complex parameter r = 
s + iX, has for each fixed value of s and A a limited set null geodesies that reach real 
3+. However, for a ribbon in the complex r-plane (i.e., a region topologically M x 7, 
with s e M and A G / = [A,„in, A,„ax]), there will be many null geodesies intersecting 
3"*". Such null geodesies were referred to as real" geodesies . More specifically, for 
a fixed s, there is a limited range of A such that all the real null geodesies intersect 
3^ in a full cut, leading to a one-parameter family of real (distorted sphere) slicings. 
The ribbon is the generalization of the real world-line and the slicings are the ana- 
logues of the spherical slicings. When the ribbon shrinks to a line it degenerates to 
the real case. We can consider the ribbon as a generalized world-line and the "real" 
null geodesies from constant s portion of the ribbon as a generalized light-cone. 

• For the case of asymptotically flat space-times, the real light-cones from interior 
points are replaced by the virtual light-cones generated by the asymptotically shear- 
free NGCs. These cones emanate from a complex virtual world-line z° = ^°(t) in the 
associated H-space. As in the case of complex Minkowski space, there is a ribbon in 
the r-plane where the "real" null geodesies emanate from. The "real" null geodesies 
coming from a cross-section of the strip at fixed s (as in the complex Minkowski 
case), intersect 3^ in a cut; the collection of cuts yielding a one-parameter family 
of cuts. The situation is exactly the same as in the complex Minkowski space ease 
except that the spherical harmonic decomposition of these cuts is in general more 
complicated. 
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Example: The (charged) Kerr metric 

Considering the Kerr or the charged Kerr metrics (or even more generally any asymptot- 
ically flat stationary metric), we have immediately that the Bondi shear tr'^ vanishes and 
hence the associated "H-space is complex Minkowski space [30, 31]. From the stationarity 
and a real origin shift and rotation, the complex world-line can be put into the form 

r(T) = (T,0,0,za), (3.12) 

with a being the Kerr parameter. The complex cut function is then 

u 

so that the angle fields are 

L 
L 
L 

Using T — s + iX in Eq.(3.13), the reality condition u = u on the cut function is that 

A = A(.,C,C) = ^<3(C,C), 

so that on the r-ribbon, A ranges between ± \/2 and the real slices from the ribbon becomes 
simply u = s/\/2. □ 

Though we are certainly not making the claim that one can in reality "observe" these 
complex worle-lines that arise from (asymptotically) shear-free congruences, we neverthe- 
less claim that they can be observed in a different sense. In the following section we 
will show that there are simple physical measurements that do determine these compex 
world-lines. 

4 Applications 

We can now explore uses of our observations concerning light-cones and their generaliza- 
tions. The first issue addressed is the application, in Minkowski space, to the Maxwell 
equations and in particular to asymptotically vanishing Maxwell fields with non-vanishing 
charge q. Specifically, we show that such solutions naturally define a complex world-line 



e{r)ia{c,c) 



(3.13) 



V2 2 

1-CC 



1 + CC' 



V2i 



l + CC 



= -V2i 



1 + CC' 



1 + CC 
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that can be identified or referred to as the complex center of charge. It is determined from 
the complex Minkowski space points where the (suitably defined) complex electromagnetic 
dipole (a combination of the electric dipole moment plus "i" magnetic dipole moment) 
vanishes. 

The analogous problem for asymptotically flat space-times (either vacuum or Einstein- 
Maxwell) is addressed with a unique world-line again arising, this time from the gravita- 
tional part with its identification as the complex center of mass. These are the H-space 
points where the complex gravitational dipole (identified as the mass dipole plus "i" an- 
gular momentum) vanishes. For the Einstein-Maxwell case there will be, in addition, a 
complex center of charge line. 

A few words of explanation in a much simpler situation might be of use. In Minkowski 
space, in a given Lorentz frame and coordinate origin, with given charge and current 
distributions (or given mass and spin distribution) , one defines the electric dipole moment 
(mass dipole) on any time slice by an space integral over the charge density (or mass 
density) times the position. By shifting the spatial origin, the dipole moment becomes a 
space-time field depending on the origin shift, 

IJ* = 1^ -qll. 

The zero values of this field determine the center of charge (or center of mass); 
By extending this idea to include the magnetic dipole moment 

'^E&zM = + i^ Al- 
and allowing the position, it, to take on complex values, we find the space dependence of 
the complex dipole moments given by 

^B&M = ^B&M - ql^c- (4.1) 
so that the complex center of charge is given by = or 

= q'^DsScM- (4.2) 

The difficulty with this construction is that it is not Lorentz invariant: the transfor- 
mations of the dipoles from one Lorentz frame to another is non-local and one does not 
obtain (in any obvious manner) a unique center of charge/mass world-line. 

We use an alternate procedure to find the different "centers of motion." Namely, the 
complex dipoles are first identified from the asymptotic solutions with interior sources: 
they are identified from the I = 1 harmonics in the tetrad components (spin-coefficient 
components) of the asymptotic Maxwell field and the asymptotic Weyl tensor, Eqs.(2.40) 
and (2.41). (See the discussion immediately after Eq.(2.43).) These quantities depend on 
the choice of the tetrad vectors at 3+. If we choose the tetrad so that the null vector 
I = l*Q determines a shear-free (or asymptotically shear-free) null geodesic congruence that 
focuses on points in complex Minkowski space (or the GR case, on points in the virtual %- 
space), we see that the associated dipole is a function (three complex components) on the 
complex Minkowski space (or H-space). The vanishing set of this function (generically) 



19 



determines the complex world-line that is referred to as the complex center of charge or 
mass. The idea is then to express the moments in terms of the complex world-line - or as 
an alternative, find the complex world-line in terms of the complex dipole. To impliment 
this (in principle straightforward) procedure is in practice rather involved, requiring se- 
vere approximations and Clcbsch-Gordon expansions of spherical harmonic products. We 
illustrate the procedure in detail with the Maxwell field in flat space and then report the 
results (obtained earlier) for the Einstein and Einstein-Maxwell cases with a minimum of 
detail. 

4.1 Maxwell Fields in Minkowski space 

Beginning with a complex world-line in Mc, z"' = £,°'{t), its family of cuts of is, as 
discussed earlier, 



with L the angle field of its null normals. Note that c has been explicitly reintroduced 
so that the cut function u, with Uret and r, have the dimensions of time. This has the 
annoying affect of causing the frequent appearance of c. 



Derivatives with respect to Uret are denoted by a prime: du^^^F = F'. 

We now illustrate how an asymptotically flat Maxwell field with non-vanishing charge 
determines a unique complex center of charge world- line, ^"(r). 
We have, first, the asymptotic solution 



(4.3) 



RemEirk 4.1 To avoid a plethora of terms involving \f2 we switch from the Bondi time 
u = ub, to the retarded time, Uret = V^ub so that 





(4.5) 



with the spherical harmonic decomposition 



4 
0' 



q + cj^lY^, + 0;,^.r4 + 



(4.6) 
(4.7) 



(4.8) 
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and physical identifications 

</.0 = 2qT^\Uret)Y,\+C-'Q'^'{Uret)Y2\^+.- (4.9) 

= -2qcrS'''{uret)Y;:'~y~'Qf'{ur,t)Y-^} + ... 

The quantities qi]^ = D^^^m — + *-^m ^^'^ Qc respectively the complex (electric 
and magnetic) dipole and complex quadrupole. 
Under the null tetrad rotation, Eq.(2.21) 

Z"^ ^ l*"" ^l-^ ^c-m" -c-m" + 0{r~'^), (4.10) 
r r 

m*'' = m^-c-n", (4.11) 
r 

n*" = n", (4.12) 

the leading Maxwell field terms transform as 

05° - 0°-2cL</);+c^L2^o, (4.13) 
0*1° = 0?-cL0°, 

The procedure to determine ^"(r) is the following: 
In the first equation of Eq.(4.13), written as 

0° = 0S° + 2cL4>l - c^Lm (4.14) 

replace the u^et (appearing in 0°, 0° and 0°) by Uret = t — ^c~^^*Fj° from (4.3), and 
for fixed r, assume that the / — 1 terms in (/)g° vanish. 

Formally, by extracting the remaining I = 1 terms in Eq.(4.14) via the integral at 
constant r, 

<L (l^lY^^dS = / (2cL0° - c2L20°)r-idS', (4.15) 

we have the exact functional relationship between the dipole gry* and the world-line ^°(r). 

Unfortunately, it is extremely difficult to get explicit relations from Eq.(4.15) and 
approximations applied to Eq.(4.14) must be used. Our basic approximation is to consider 
the £,°'{t) to be of the form ^°(r) = (r, ^'(r)) (using a r re-scaling of the form r — ^ F{t)) 
with both and the rf to be "small." We retain only terms up to second order and 
harmonic expansions up to Z = 2. 

Writing out Eq.(4.14), 

0° = 0S° + 2cL0° - c^L^^o, 
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using Eqs.(4.9), with Uret = t — and L = c ^C{''')Y-^i, by omitting cubic 

terms including , then using the first two terms of the Taylor series 

F (^r - ^c-if ^ = F(r) - ^c-^f r°f^'(r), (4.16) 

and the Clebsch- Gordon expansions of the products of the spherical harmonics (see Ap- 
pendix A) we finally have (after simplification) for just the / = 1 harmonic terms 

qv" = qe - ^\ec-W'eak + ^c^^g^f • (4-17) 
The first thing we notice is the linear relation: 

V'{r)^e{r). (4.18) 

This can be fed back into Eq.(4.17) in either of two ways resulting in either of the 
relations: 



■q = ^ - i^c ^ ^ tuk + -^q c Qc i , (4.19) 

s — '? +*2 — 10^^ vc s J 

which determine the complex dipole in terms of the complex world-line or the world-line 
in terms of the complex dipole. 

Note that though all the expressions are functions of r, the r can be replaced by 
Uret with no other changes needed due to our approximation scheme. 

4.2 Asymptotically Flat Space-Times 

Turning now to the Einstein (or Einstein-Maxwell) case, we basically repeat the procedure 
used in the Minkowski space Maxwell field example. 

We begin with an unknown complex world-line in "H-space, z"" — C(t), to be deter- 
mined by the existing asymptotically fiat space-time. Its family of cuts and null normal 
angle field of is, as discussed earlier (c.f., (2.36), (2.33), etc.), given by 

u = c-1G(t, C,C) (4.20) 

= C-lr(r)fa(C,C) + C-^i?i^2(r(r),C,C)=C-lr(T)/a(C,C)+C-^i/i^2(T,C,C), 

LiuXX) = c-ig.G(r,C,C)I..T(u,c,C)=c"'^'W^H-6c-'^''W^H.+--- (4-21) 
a"(r,C,C) = g^,)G(T,C,C)-24f^"(r)y4.-f ... (4.22) 

We now show how a given asymptotically fiat space-time determines the complex center 
of charge world- line (r) . 
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Returning to the "peeling" theorem: 





= i>y 


+ 0{r- 


-% 


Ipl 




-4 + 0(r- 


-% 




= 


-^ + 0{r- 






= 


-^ + Oir- 




tp4 









with the transformation law of the leading terms under a null rotation, 



= v;i-4cLVjf + 6c2l2^;°-4c3l3^*° + c4lV4, (4.23) 

= i/,0-3cLV2+3c2iV3 -c^i^V-", (4-24) 

V2*° = V2-2cLV3+c'iV2, (4.25) 

= ^^l-cLi,l (4.26) 

= i>l (4.27) 



we can then determine the transformation law for the physical quantities that are identified 
in the following harmonic components. 

The (truncated) harmonic expansions with their (approximate) physical identifications 

are 

< = <''yiii+- (4.28) 

^? = VWi + - 

^0 = VI/ _ g2-o - c-VO(^O)' 
= f = ^fO + ^r'yo + ... 
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and 



i/jg'"' — approximately, the quadrupole (4.29) 



^C(g™.) = Dinars) J' ^-^^l' (4-30) 

^1 = -^^iDUass)+^c-'J') (4.31) 

* = V^° + gV + c-iaV° (4.32) 

^ = ■* = «'" + 'J''y/J + ... (4.33) 

2\/2C 6G* 

= -Mb—^ tP'Yi, + ... Mass Aspect 

c 

Ms = Bondi mass (4.34) 

2v 2G 

P' = Bondi linear momentum (4.35) 

6G 

cr" ^ 24f ^^2^ + • • • Bondi asymptotic shear. (4.36) 

= (e]^■+^^}^■) = _|_(g;^'4^ + ^Q^;J (4.37) 

Remark 4.2 T/ie relationship between ipi^ and the mass dipole, -D^^^^^j, and angular mo- 
mentum, J' , is usually considered to be more complicated than stated here, often involving 
quadratic terms in the Bondi shear [32]. The trouble is that there are disagreements in 
these quadratic terms in the different versions. We have simply left them out here and 
note that in our approximations the disputed terms do not appear. 

Concentrating on the transformation of the dipole, Eqs.(4.23)-(4.27), our focus is only 
on Eq.(4.24), which can be rewritten as 

^l,^^ = + 3cLi/'2 - 3c2lV3 + c^^VS- (4-38) 
We first note that the I = 1 term in is proportional to the complex gravitational 
dipole, -C)c(gr ai.)(")- procedure is now to replace all the us that appear in Eq.(4.38) 

by uret = {V2c)-^^''{T)la{C,0 + (^) "'^f;^2 (^ (t) , C, C) (1-6., Eq.(4.20)), and remember 
that L also depends on (r) . Then, from our basic assumption, we take the / = 1 term in 
ipl'^ to vanish and finally extract the I — 1 harmonic coefficients from Eq.(4.20). Formally, 
this is done via the integral expression 

(f -iplY^^^dS ^ <L {3cLij°-3c^L^^^ + c^L^tpl)Y{-^dS, (4.39) 

which on the left-side contains r and the dipole ^c(grat)) while the right-side contains 
T, the unknown world- line S,°'{t) and the Bondi shear, cr". 

Though in principle this equation should allow us to establish the relationship between 
^C{grav) ''^^'^ i''') ' practice this is not possible: we must to return to Eq.(4.38) and use 
harmonic and Clebsch- Gordon expansions with severe approximations and finally collect 
the I = 1 terms directly. 



V2 
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4.2.1 A Poor Approximation: Results 



We first describe a prelirainary procedure for extracting tlie compiex world-fine from an 
asymptoticaffy flat space-time rather explicitly. The foflowing approximations are used: 
the Bondi mass is taken as zeroth order while all other variables are first-order with the 
calculations done keeping terms up to second order. In this preliminary version, the 
harmonic expansions keep only the I = (0, 1) terms. This implies the severe condition that 
the Bondi shear (cr'^) be taken to be zero. Later this condition is relaxed. 
Via these approximations Eq.(4.38) becomes 

^PUu, CO = + 3cL{t, CO^iu, C, C). (4.40) 

Using the retarded time u, et — V^u instead of the Bondi time, then replacing all the 
UretS by Uret — T — "^c^ (t) Yj^^j and finally Taylor expanding with (4.16), Eq. (4.40) 
becomes (with physical identifications inserted): 



r,\r)Y,\ - ^c-Vr'(r)eyi°n\ = -^-^Mne{r)Y,\ - 1^ ir)Y,'^,) . (4.41) 



Finally after Clebsch-Gordon expansions and the use of the (linearized) Bianchi iden- 
tity, Eq.(2.44), 

V2tlj°' = 2c-f'Y,\ ^ = V2c^'* = -^^^P\ (4.42) 

the three complex I = 1 coefficients of (4.41), (with ^'^ = + yield 

(D^mass) + ^C-'J') = MBi^R + ^^i) " e^.^c-' P^l^,^ - ^j) , (4.43) 
or the pair of real equations, 

DUass)ir) = MBeB.{r)+c-h,,kP'^j, (4.44) 

j'iT) = cMB^IiT)+e,,kP'(Ur)- (4.45) 

From Eq.(4.42), we immediately get the kinematic definition of the Bondi momentum 
in terms of the complex world-line. In addition, we have the conservation of angular- 
momentum which arises from the reality of the mass aspect, ^, Eq.(4.33): 

P' = Dl'^ass)=MBeR +c''e,jkMB{CReiy, (4.46) 
j'' ' = 0. (4.47) 

There are several things of significance that should be pointed out here. 

• If the higher gravitational moments and electromagnetic terms were included, these 
results, Eqs. (4.44)- (4.47), would all be augmented by further terms. In particular 
there would be a non-vanishing angular-momentum flux. See below. 
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• In the expression for the angular momentum there are two terms, the second be- 
ing the conventional orbital angular momentum while the first has been identified, 
via the Kerr metric and the charged Kerr metric [33], as the intrinsic spin angular 
momentum. 

• The mass dipole contains the conventional Af plus a momentum-spin interaction 
term that creates a spin- velocity coupling contribution to the linear momentum. 

These results - basically kinematic, aside from the conservation of angular-momentum 
- have been derived, with severe approximations, by associating the idea of a complex 
center of mass curve with a complex curve in ?^-space. 

In the same vein (with the same approximations), we obtain the dynamic law for the 
motion of the real part of ^'^ (i.e., ). From the second Bianchi identity, Eq.(2.45) 

V2^0' = -cStAO ccrV2 ^ i^2=0 (4.48) 
we have that Mb and are constant, i.e. conservation of energy and momentum: 

M's = P" = 0. (4.49) 

4.2.2 A Better Approximation (with Bondi shear): Results 

For a more accurate description/determination of the complex world- line associated with 
a given asymptotically flat Einstein (or in the following subsection, Einstein-Maxwell) 
space-time we restore, in the calculations, the Bondi shear and include the effects of the 
Einstein-Maxwell equations. Rather than redoing the calculations from the beginning, 
using the same procedures as in the previous section, we simply give the final results. 
The approximations are basically the same: the Bondi mass is zero-order, while all other 
variables are first -order; in the calculations only quadratic terms are retained. However 
the harmonic expansions now include the I = (0,1,2) harmonics. In addition, since a 
Maxwell field (with non- vanishing total charge) is allowed, we have, not only the complex 
center of mass line, z°- = C(t), but at well the complex center of charge line. It is denoted 
by 2° — ?7°(t). In general the two line are different, though in special circumstance they 
can coincide. 

The idea is to start with Eq.(4.40), use the known expression for ilji{uretX^C) ^^nd 
^{ureti C: C)i i^i tcrms of the physical gravitational moments, then replace every Uret by 

uret = ^c-iX(r(r),C,C), (4.50) 

= r - ^c-ir(r)r° + V2c-ir^-(r)y°,. + ... 

set the I — 1 coefficients of tAi" to zero and finally extract the I = 1 coefficients from 
the entire equation - a long process involving repeated Clebsch- Gordon expansions (c.f., 
[15, 7]). This process leads to 
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<^ = -6V2Gc-2(i?|„_)+zc-ij') (4.51) 



4G If 3 -rrrilff 



5c^ ^Grav VGrav'^U* 



or from the real and imaginary parts, 



4G 



ass 



5^2 WMassSK + ^/SpinS/J g^g "^Mass Vspin^'j* 
4G 3 

the definition of the mass dipole and angular momentum in terms of the complex world- line 
The kinematic definition of the linear momentum and the angular momentum conser- 
vation law are then found by extracting the I = 1 harmonics from the Bianchi identity, 
Eq.(2.44), (the evolution equation for V'l) 



^0 = -g?/,o + 2ctV, 



0„;.0 
3' 



V'?' = -^cg* + ^c5^a + 3a°g(a'). (4.52) 



leading to 



Then inserting the expressions for Dl sand J' we obtain, from the real part, an 



expression for the linear momentum 



— D\mass) (QMa'ssQs'pin ) ^J^i^ (4.53) 



12G 

(mass) l^^Mass^Spin^ 

= A/bCh+^P'. (4-54) 
3G 



= c-l(e|P^)'6,H--.(P^Qr.s)'-^(0MlsCfl + Qf;i„el)' 



^6 (^'^MassQs'pinj ^i^j' 

and, from the imaginary part, the angular momentum conservation law: 



J*' = (Flux)\ (4.55) 
4G 3 

r = MBcC, + P^eRel,^-^P'Q%L~-{QXsseI-Q%.JR). (4.56) 

(Flux)^ = ^{q^^^L + QI!^LQ^X)'^^>'- (4-57) 
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Finally, from the Z = 0, 1 parts of the Bianchi identity, Eq.(2.45), the evolution equation 
for the mass aspect 

or 

VP' ^ ^aO'^O', (4.58) 
c 

we obtain both the energy loss expression and the evolution of the momentum, i.e., the 
equations of motion. 

The energy (mass) loss equation, from the I — part, is 

G 



M' — - — ( O*-' " O*-' O*-' " D'^ "' 

iVJ^B — ^^j \^'VMass'VMass + ^SpinVSpin 



the known quadrupole expression, while the momentum loss equation, from the I = 1 part 
of (4.58) becomes a version of Newton's second law: 

P" = FLoii: (4.59) 

-^recoil — -j^g^e ^^Spin^Mass ^Mass^Spin/*^2'fc ' 

Finally substituting the P' from Eq.(4.53), we have Newton's 2nd law of motion; 

MbCr = F^,,, - Mj^a - = (4-60) 
4.2.3 Results for Einstein-Maxwell Space-times 

The calculations that were performed earlier for the vacuum GR case can be extended 
to the Einstein-Maxwell case with considerably more effort. Rather than going into the 
details we will simply present the main results. The Maxwell field considered has only 
charge and dipole terms: with a bit of effort quadrupole terms could be included. The 
main change needed is the modification of the asymptotic Bianchi identities to include the 
Maxwell field: 

40 = _g^o + ^o^o^fc^o^o^ (4.61) 
V;? = -gi/.^ + 2crV!^ + 2/c0?,^°, (4.62) 
k = 2Gc~^. (4.63) 

where the fields (t)\ and (f)^ are given by Eqs.(4.9). 
The result of the calculations are 

= -^-^{DUass)+^c-'r), (4.64) 

{dUss)+^c-'j') - M^e + '-^rp' -^^,^rn.k^t'' ^^.p'qZl (4.65) 

15^6 " *^Giav g^2^ '^Grav * ^'"j'^'^Grav^^Grav 
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DUass) = MB^},-c-'P^^'}e,,.-^{enQZss + ejQZ'in) (4-66) 



and 

= cMBei + enP'e,,,-^^[QZssei-QZLeH) (4.67) 



From the Bianchi identity, Eq.(4.62), we have 

-'^(mass) + *C J - -f" + + ^^(^GravyGrav^jfe* ^ 3^5^^ ^ ^ 

or, from the real part, 

P' = MBCk-%rtl(-c-\P=ei)'euj. + ^,{eiv'i^ -e^vT)' (4.68) 

5^2 \^?i?/^Mass S/^^Spin^ ^6 y^^Spin'^MassJ ^J^*' 

and imaginary part 

r' = (Flux)' (4.69) 

\fmyi} — [VbVr J ^jkz + ^.^g l^yMassVMass + VSpinVSpinJ <^fcj* Sc^^' 

Alternatively, we can consider the term ffr??/" as a contribution to the total angular 
momentum rather than the flux. Using this alternative definition of angular momentum, 

2(72 

-Jt - J + ^Vi , 
where J* is given by Eq.(4.67), leads to a modified flux law: 

= (Flux)^ (4.70) 
(Flnx)^ = ^ + ^fr^r) + ^ {qZssQ'^L + Q^tinQr/n) ^>^3^- 
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Finally from the Bianchi identity, Eq.(4.61), we have the mass loss and momentum loss 
equations: 

S (^QMass^Mass + Qspin^Spin j " 1^ [Vr Vr + Vi Vi ) , (4.71) 

" \^ (Qspm^Mass " QMass^SpinJ ^Ikr + ^ (?7/ ^/fl - Vr V'l ) <^jH- (4.72) 

4.3 Interpretations 

There are a variety of comments to be made about the physical content contained in the 
above kinematic and dynamical relations: 

• A subtle (but not essential) comment for completeness should be made. We have 
used the symbol gry* for the complex electromagnetic dipole and 77" for the complex 
center of charge. The two rjs are related by a non-linear term, given explicitly by 
(4.19)). Rather than go into a detailed explanation of our usage which could be 
confusing, we note that the effect is that a small quadratic term in is missing. 

• The first term of Eq. (4.68), is the standard Newtonian kinematic expression for 
the linear momentum, MbS,%' ■ This is followed by a spin-momentum coupling term 
of the form x ^)'. 

• The further term, —^c~^q^rf^\ which is a contribution to the linear momentum from 
the second derivative of the electric dipole moment, gry]^, plays a special role for the 
case when the complex center of mass coincides with the complex center of charge, 
rf =£,"'■ In this case, the second term is exactly the contribution to the momentum 
that yields the classical radiation reaction force of classical electrodynamics, [34] 

^c"3<?er- (4.73) 

In this special case we have a rather attractive identification: since now the magnetic 
dipole moment is given by D\[ — q£^} and the spin by 5' = MbcS^}, we have that the 
gyromagnetic ratio is 

^ _\S^ ^ Mbc 
^ \DU q ' 
leading to the Dirac value of the g-factor, i.e., g — 2. 

• Many of the remaining terms in P*, though apparently second order, are really of 
higher order when the dynamics are considered. Others involve quadrupole interac- 
tions, which contain high powers of c~^. 

• The complex electromagnetic dipole moment, given in general by 
becomes, when the world-lines coincide, D^^^^^j — qi^R + «^1). 
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In the expression for J* (4.67) we have already identified, in the earlier discus- 
sion, the first two terms, S'^ — MbcS,^ and Ms^.'^^'CR^ikj as the intrinsic spin an- 
gular momentum and the orbital angular momentum (i.e., x i^) respectively. 
An interesting contribution to the total angular momentum comes from the term, 
^c^^q^riY — ^c~^qD'^j, i.e., a contribution to the total angular momentum from a 
time- varying magnetic dipole. A question arises: is this an observable prediction? 

Our identification of J* as the total angular momentum in the absence of a Maxwell 
field agrees with most other identifications (assuming our approximations) [32]. Very 
strong support of this view, with the Maxwell terms added in, comes from the flux 
law. In Equation (4.70) we see that there are four flux terms (more arise if we in- 
cluded electromagnetic quadrupole radiation): the first and second come from the 
Maxwell dipole flux, while the third and fourth are the gravitational quadrupole flux 
terms. The Maxwell dipole part of the flux is identical to that derived from pure 
Maxwell theory [34]. We emphasize that this angular momentum flux law has little 
to do directly with the chosen definition of angular momentum. The imaginary part 
of the Bianchi identity. Equation (2.45), is the conservation law. How to identify the 
different terms (i.e., identifying the time derivative of the angular momentum and 
the flux terms) comes from different arguments. The identification of the Maxwell 
contribution to total angular momentum and the flux contain certain arbitrary as- 
signments: some terms on the left-hand side of the equation, i.e., terms with a time 
derivative, could have been moved onto the right-hand side and been called "flux" 
terms. Our assignments were governed by the question of what terms appeared most 
naturally to be explicit time derivatives (thereby being assigned to the time deriva- 
tive of the angular momentum) , or which terms appeared to be physically more likely 
to be an angular momentum term. 

The angular momentum conservation law can be considered as the evolution equation 
for the imaginary part of the complex world line, i.e., ^}(ui.ct). The evolution for the 
real part is found from the Bondi energy-momentum loss equation. 

The Bondi mass, Mb = — 2^20 ^ ' ^^'^ original mass of the Reissner-Nordstrom 

(Schwarzschild) unperturbed metric, Mrn — "5^20^°" (i-^-' ^ ~ ^ harmonic of 

1P2) differ by a quadratic term in the shear, the I — part of a'^a'^ . This suggests 
that the observed mass of an object is partially determined by its time-dependent 
quadrupole moment - if it exists. 

In the discussion of the Bondi energy loss theorem (4.71), we saw that we can relate 
^'■^ (i.e., the I = 2 shear term) to the gravitational quadrupole by 

to obtain the standard quadrupole energy loss. 

The Bondi mass loss theorem with electromagnetic dipole and quadrupole radiation 



31 



becomes 



M'b - -fj(Q^::^L + mZ^X) (4-75) 

with the first term the conventional gravitational radiation with the second and third 
terms the electromagnetic dipole and quadrupole radiation loss. The momentum loss 
with electromagnetic quadrupole contributions becomes: 

pj7 _ ( ri^jiii r\^jiii ri^jiii r^iim \ i 9^ ( kn jii kii jii\ 

~ 15c6 V^Spin^Mass '^Mass'^Spin J ^Iki + yll Vr Vr Vi J ^jki 

I ^ ( r^ljiii r)kjiii _ r^ljiii r)kjiii\ 

• There are several things to observe and comment on concerning Eqs. (4.60) and (4.59): 
Returning to the case when the complex center of mass and center of charge coincide 
the resulting equations of motion for the world-line, namely: 

we observe and stress that, aside from several extra terms, these equations coincide 
with the Lorentz-Dirac equations of motion - this includes the hyper-acceleration 
term and the mass loss term, M^'^]^. This result follows directly from the Einstein- 
Maxwell equations. There was no model building other than requiring that the two 
complex world lines coincide - a strong condition, perhaps equivalent to a point par- 
ticle assumption. Furthermore, there was no mass renormalization; the mass was 
simply the conventional Bondi mass as seen at infinity. Further structures (e.g., 
spin) could remain. The problem of the runaway solutions, though not solved here, 
is converted to the stability of the Einstein-Maxwell equations with the "coinciding" 
condition on the two world lines. If the two world lines do not coincide, i.e., the 
Maxwell world-line is formed by independent data, then there is no problem of un- 
stable behavior. This suggests a resolution to the problem of the unstable solutions: 
one should treat the source as a structured object, not a point, and centers of mass 
and charge as independent quantities. Alternatively, it might be possible that the 
extra terms in the equations might stabilize the equations. It is however hard to see 
how this could be demonstrated. 

• The i^rccoii is the recoil force from momentum radiation, other force terms could be 
considered as gravitational radiation reaction. 

• There are alternative perturbations schemes that use variations of the procedures 
used here. An example is the determination of the gravitational radiation in a 
Schwarzschild or Reissner- Nordstrom space-time induced by a time-dependent Maxwell 
dipole radiation field. The physical identifications agree with those found in the 
present work. For instance, the perturbations induced by a Coulomb charge and 
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general electromagnetic dipole Maxwell field in a Schwarzschild background lead to 
energy, momentum, and angular momentum flux relations [35]: 



M'r: 



B 







(4.76) 



P' 




J' 



-ki 



— ( " '4- " ' ^ 



^ijk^ 



all of which agree exactly with predictions from classical field theory [34] . 

The familiarity of these results and those in the full text act as an exhibit in favor of the 
physical identification methods described in this work. That is, they act as a confirmation 
of the consistency of the identification scheme. 

5 Discussion &; Conclusion 

We have studied and described a variety of geometric structures, all occurring within the 
confines of classical special and general relativity, which strongly resemble ideas in other 
areas of physics. Our complex-conjugate method for describing a real, twisting shear- 
free, or asymptotically shear-free NGC places the congruence's caustic set (in Minkowski 
space, interpreted as its source) on a closed curve propagating in real time; or in more 
suggestive parlance, a classical string or world-tube. The dual description, or holomorphic 
method, constructs a complex NGC (a complex light-cone congruence) whose apex is on 
complex world-line in either complex Minkowski space (the shear-free case) or "H-space 
(the asymptotically shear-free case). When one imposes a reality structure on the null 
geodesies of this congruence (i.e., asks that they intersect the real asymptotic boundary), 
the source becomes a real two-dimension (complex 1-dimensional) open world-sheet, also 
oddly reminiscent of string theory. 

Furthermore, the role of 3^ as the object which interpolates between these two dual 
descriptions is reminiscent of the conjectured holographic principle [2, 3]. In particular, 
we can think of a real, twisting asymptotically shear-free NGC as determining data on 3^7 
(i.e., the complex-conjugate construction), which in turn acts as a "lens" into 'H-space, 
which serves as a virtual image space for the real space-time, where the real twisting 
NGC becomes a complex, twist- free NGC. Hence, one is tempted to refer to 3^ as the 
holographic screen for some application of the holographic principle to classical general 
relativity. This should be contrasted against the most famous application of the holo- 
graphic principle: the AdS/CFT correspondence [4, 5]. Here, the AdS boundary acts as a 
holographic screen interpolating between type HB string theory in AdS^ x and A/" = 4 
super- Yang-Mills theory in Minkowski space-time. It is interesting that we have found 
structures so closely resembling the underlying holographic principle, but which involve 
nothing like extra dimensions or supersymmetry. 

Our observations also raise a series of potentially interesting questions related to flesh- 
ing out the alluded-to connections between general relativity and more ambitious theories. 
For instance, is it possible to write down a sigma model for the embedding of the open 
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world-sheet (the holoniorphic method) into 'H-space which yields the real cuts of 3+? Also, 
can the role of the future asymptotic boundary as a holographic screen be made any more 
precise? Furthermore, we have seen (with several examples) how the geometry of the 
virtual image space (?^-space) allows us to make physical identifications from data on the 
holographic screen. This indicates that our virtual image space represents the "physical 
information" side of some holographic principle, in the same way as the conformal field 
theory (CFT) side of the AdS/CFT correspondence. 

It should be noted that in 't Hooft's original work connecting gauge theory and string 
theory in the planar limit, no supersymmetry was introduced [1]; an extra dimension for 
string propagation does enter for anomaly cancellation, though. It could be possible that 
the analytic continuation of to the holographic screen 3^7 in our investigation serves an 
analogous purpose, adding the degrees of freedom necessary to solve the good cut equation 
and construct the virtual light cone congruence. However, it is unclear whether further 
degrees of freedom (in the form of additional dimensions) would be needed if our notion 
of the source of "real " null geodesies as an open string in H-space were taken seriously 
(i.e., if one attempted to quantize the theory of such "strings"). 

Other (perhaps less abstract) questions immediately arise. Since we obtain many of 
the standard classical mechanics kinematic and dynamic equations, is it possible that the 
standard quantization could shed light on the difficult issues of quantum gravity? Since 
the material described here is closely related to Penrose's twistor theory (via different 
versions of the Kerr theorem relating shear-freeness and twistors; c.f., [36, 7, 31]) is there a 
twistorial version of the present findings? Though the material presented here deals with 
the apparent or "virtual" motion of compact sources as viewed from a great distance, can 
it be modified or generalized to treat interacting bodies? Can or should these complex 
world-lines with dynamic evolution and spin structure be taken at all seriously - and in 
what context? 
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A Tensor ial Spin-s Spherical Harmonics 

Throughout the text, we have made use of tensorial spin-s spherical harmonics, Yf'^ ■ (C, (^), 
to expand analytic functions on the 2-sphere. Although these are closely related to the 
usual spherical harmonics familiar from other arenas of physics, the incorporation of spin 
may be new to some readers. This appendix reviews the basic properties of the spin-s 
harmonics and includes some of the Clebsch-Gordon expansions of their products used in 
our calculations. For the original treatment of this material, see [24]. 

Before introducing these functions though, we recall the definition of the 3-operator 
(also used throughout the text), which acts as a spin- weighted covariant derivative on 5^ 
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[27]. Suppose f : C has spin- weight s; then 



dC, dC 

where P is taken to be any conformal factor on the 2-sphere (in the text, we always 
take P = 1 + (^). We see that an apphcation of the S-operator raises spin- weight by 1, 
while 8 decreases it by one. As we will see, this property allows us to move among the 
spin-weighted tensorial spherical harmonics simply by applying differential operators. 

A.l Spherical Harmonics 

Some time ago, the generalization of ordinary spherical harmonics YimiCTC) to spin- 
weighted functions (s)Yim{CX) (^-g-, [37, 27]) was developed to allow for harmonic ex- 
pansions of spin-weighted functions on the sphere. We have used instead the tensorial 
form of these spin-weighted harmonics, the so-called tensorial spin-s spherical harmonics, 
which are formed by taking linear combinations of the (s)Yim (Ci C) [24]: 

i...k — 2^ i...k (s)'^lm, 

where the indicies obey ]s] < I, and the number of spatial indicies (i.e., i...k) is equal to 
I. Explicitly, these new spin- weighted harmonics can be constructed from first defining 

i...k- 

We impose our special 2-sphere Bondi null tetrad: 

= -^(P, C + C" -«(C-C), -1 + CC), (A.2) 
^ (p, ~(C + C), *(C-C), 1-CC), 



\/2P 

= ^(0, l-C^ -*(i + C'), 2C), 
p = 1 + CC 

We then project these to their covariant duals, and take the spatial parts to obtain the 
1-form components, along with a useful additional form Ci = li — nf. 

h = -^(C + C -«(C-C), -1 + CC), (A.3) 
^ (C + C, -«(C + C), -1 + CC), 



rrii 



V2P 

^ (l-c^ -*(i+c^), 2c) 



V2P 

From this, we define ^ j, as [24] 

Yi\...k = rnimj...mk, (A.4) 
YfL.k = rhifhj...fhk. 
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The other harmonics are determined by the action of the S-operator on the harmonics 
we have defined in (A. 4). In particular, it can be shown that 

Yi\..k = {Yi\...k) , (A.5) 



i...k ^ " Vl I.. 



k) 



We now present a table of the tensorial spherical harmonics up to Z = 2, where we 
truncate all of our expansions in this paper. Higher harmonics can be found in [24]. 



• / = 

• / = 1 



TO,; 



/ = 2 



Yu^ = ™, 



Y2ij = m^rrij Fjij = " (cjWj + rriiCj) Y^^^ = 3c,Cj - 25ij 
Y2if = fhirhj Y^ll = - {ciffij + fhiCj) 



In addition, it is useful to give the explicit relations between these different harmonics 
in terms of the 5-operator and its conjugate. Indeed, we can see generally that applying 
S once raises the spin index by one, and applying 3 lowers the index by one. This in turn 
means that 

^Yi\,„k = 0, 

^YrLk = 0. 

Other relations for Z < 2 are given by: 

sro = ~2Y,\ 
8F° = -2y-i 



ax^2 _ T^l 

S2y2 _ yO 
•J ''2ij — ^2ij 

= -^Y,% 

Finally, due to the non-linearity of the theory, throughout this review we have been 
forced to consider products of the tensorial spin-s spherical harmonics while expanding 
non-linear expressions. These products can be expanded as a linear combination of indi- 
vidual harmonics using Clebsch-Gordon expansions. The explicit expansions for products 
of harmonics with / = 1 or / = 2 are given below (we omitt higher products due to the 
complexity of the expansion expressions) . Further products can be found in [24] . 
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A. 2 Clebsch-Gordon Expansions 

• ^ = 1 with / = 1 

• / = 1 with / 2 

^li^2ij — ^3ijk 

Yi'^2% = lY^^S,k - lY,)S,k - Iy,\S,, + -J= {e,kiY,]i + e,,iY^^{) + '^Y^^^ 

Y^^Y^k - ^Y,%k + ^Y,°,S.k - \Ylk, + ^ {e,kiY,l + e^uY^) - ^Y,%, 
Y^^Y^.k = -lYl,5,k + \Yl^5.k + \Ylk5^l ~ ^ {e.uY^.i + e^jiY^ki) + ^Y^.^k 

vl v-l _ 3 4- a "'^ vO a yO J- c ^ "'^ V"0 

^2ij-'lA; - Y^^li^i*; + Y^^lj^ifc ^ g^lfe"ii \^]kll2il + ^ikll2l]) ~ 3Q^3ij7c 

'-2ij'-lk — " V^2ij J^ljt j 
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